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Abstract 

Pressure data measured on the flexible OAST 
ARW-2 wing are compared with results calculated 
using the transonic small perturbation code 
XTRAN3S. A brief description of the analysis Is 
given and a recently-developed grid coordinate 
transformation Is described. Calculations are 
presented for the rigid and flexible wing for 
Mach numbers from 0.60 to 0.90 and dynamic 
pressures from 0 to 1000 psf. Calculated and 
measured static pressures and wing deflections 
are compared, and calculated static aeroelastlc 
trends are given. Attempts to calculate the 
transonic instability boundary of the wing are 
described. 


Nomenclature 

c airfoil chord. In 

Cl wing lift coefficient 

Cm wing pitching moment coefficient 

about the wing root leading edge 
C D pressure coefficient 

C p critical pressure coefficient 

cf. wing reference chord, in 

DN value at downstream boundary of finite 

difference grid 
LE leading edge value 

M freestream Mach number 

OB value at far spanwise boundary of 

finite difference grid 
q dynamic pressure, psf 

r function defining Instantaneous 

position of wing surface, z = r(x,y,t) 
t time, sec 

TE trailing edge value 

TIP wing tip-value 

UP value at upstream boundary of finite 

difference grid 

V free stream velocity, In/sec 

x, y, z coordinates of a right hand Cartesian 
system with origin at wing root 
leading edge, positive x In downstream 
direction, y in spanwise direction, 
and z up, in 

Zt wing-tip displacement, In 

a wing root angle of attack, degrees 

y ratio of specific heats 

"n fraction of semispan 

e t wing-tip rotation, positive leading 

edge up, degrees 


5, n, c transformed coordinates in x, y, and z 
directions respectively 
p density, slugs/ft-* 

,x transformed time, t = t 

perturbation velocity potential 


Introduction 

The transonic speed range is a critical 
region for many aeroelastlc phenomena such as 
flutter and divergence. In the past, analytical 
methods have been unable to predict accurately 
the nonlinear transonic aerodynamics and 
analysts have used linear theory for estimates 
of transonic aeroelastlc behavior. Reliance was 
placed primarily on tests of scaled aeroelastic 
models and on flight tests of prototype aircraft 
for flutter clearance and aeroelastic 
deformation effects. Recently, considerable 
progress has been made in calculating steady 
transonic flows about aircraft using finite 
difference methods to obtain numerical solutions 
of the flow equations. Significant progress 
also is being made toward developing finite 
difference methods for unsteady flows which 
eventually may lead to accurate transonic 
aeroelastic analyses. 

For two-dimensional flows, methods based on 
the transonic small perturbation (TSP) equation 
have been developed and extensively applied (see 
ref. 1-2 for example). These methods have been 
extended to Include viscous effects, 3 * 4 non- 
isentropic effects, 5 and wing-canard 
configurations. 6 For three-dimensional flows, 
the XTRAN3S program has been developed 7 by the 
Boeing Company under USAF contract. It treats 
an Isolated planar wing including aeroelastlc 
deformation effects and unsteady motion. The 
XTRAN3S program has been implemented on the 
Control Data Corporation VPS 32 computer at the 
NASA Langley Research Center. A variety of 
applications are being made In order to evaluate 
Its applicability to several types of wings. 

Several organizations have expended 
significant effort to measure static and 
oscillatory pressures on wings at transonic 
speeds for use In evaluating computational 
methods and to Improve the understanding of 
unsteady transonic flows. For example, 
configurations tested at the Langley Research 
Center Include a clipped delta wing with a 6 
percent thick circular arc airfoil section, 6 an 
advanced transport wing with several oscillating 
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controls, y an oscillating rectangular wing with 
a supercritical airfoil, u and a flexible 
supercritical wing from the DAST ARW-2 
vehicle. 11 

Pressures for the rectangular wing have 
been calculated using XTRAN3S. 12 Calculations 
have also been made for the F-5 u > iy and LANN 
wings 15 which were tested at NLR, 15 « b and for 
the AGARO SMP tailplane model 1 ' tested by the 
British RAE. 1B 

These calculations give an indication of 
the capability of XTRAN3S to predict steady and 
unsteady pressures on rigid oscillating models. 
However, XTRAN3S also includes the capability of 
computing static aeroelastic deformation and 
transient time history responses for flutter 
analyses. The test of the second Aeroelastic 
Research Wing (ARW-2) in NASA's Drones for 
Aerodynamic and Structural Testing (OAST) 
program in the Langley Transonic Dynamics 
Tunnel 11 provides bn opportunity to evaluate 
this capability. Comparisons of some of the 
data from this test with aeroelastic XTRAN3S 
calculations are given in this paper. 

A brief description of the DAST ARW-2 wing 
and the wind tunnel test is first presented. 

Then overall descriptions of the XTRAN3S program 
and a recently-developed grid coordinate 
transformation are given. Calculated and 
measured results are presented for static 
pressures and for wing-tip deflections. A brief 
trend study of static aeroelastic effects is 
described, and some efforts to calculate the 
experimental transonic instability boundary are 
discussed. 


Description of Wing and Test 

The right wing panel of the DAST ARW-2 
flight vehicle was tested in the Langley 
Transonic Dynamics Tunnel (TDT) as shown in 
Figure 1. It was mounted on top of a half-body 
fuselage. The wing has a supercritical airfoil 
section, an aspect ratio of 10.3, and a 
leading-edge sweepback angle of 28.8°. The wing 
thickness is 15, 12, and 11 percent chord at 
n = 0.071, 0.426, and 1.00 (the root, trailing 
edge break, and tip stations) respectively, with 
a linear variation of thickness between these 
span stations. The wing planform and overall 
arrangement are shown in Figure 2. The wing has 
two Inboard control surfaces and an outboard 
control surface which are hydraulically 
controlled. For this test, the inboard surfaces 
were not used but were hydraulically held at 
zero deflection. The outboard surface was 
deflected both statically and dynamically in 
order to excite the model. 

The design condition for the wing Is 
M = 0.80, C[_ = 0.53, and q = 125 psf. The 
ARW-2 wing was designed to require active 
flutter control within its flight envelope. 

The primary instrumentation consisted of 
191 pressure transducers, 10 accelerometers, and 
wing root strain gages. The orifice rows of 
upper and lower surface pressure transducers 
were located at n = 0.274, 0.476, 0.599, 0.707, 
0.871, and 0.972 which correspond to the 


dimensional locations shown in Figure 2. Steady 
and unsteady pressures were measured for a large 
number of test conditions in the TDT using Freon 
as a test medium. The test conditions at which 
pressure data were obtained are shown in Figure 
3. Data were obtained at Mach numbers of 0.6, 
0.7, 0.8, 0.85, and 0.88 and at dynamic 
pressures of 100, 200, and 300 pounds per square 
foot (psf). At each tunnel condition static 
pressure data were obtained for wing root 
angles-of -attack of -2 to 4 degrees. Some of 
the higher angle-of-attack values were 
eliminated at the higher dynamic pressures due 
to wing loading limits. Unsteady pressure data 
were obtained at wing angles-of-attack of 0 and 
2 degrees for control surface oscillation 
amplitudes of 1, 2, and 3 degrees and 
frequencies of 5, 15, and 20 Hz. 

An unsual wing instability was encountered 
in the wind tunnel test. This Instability 
boundary is shown in Figure 4 and occurred at a 
nearly constant Mach number of about 0.90 for 
dynamic pressures from near 50 psf to above 300 
psf. The frequency of this Instability ranged 
from about 8.6 Hz at the low dynamic pressure to 
approximately 13 Hz at the higher dynamic 
pressure. The wing instability motion was 
dominated by the wing first bending mode which 
had a measured frequency of 8.3 Hz at wind-off 
conditions. The instability was found to be 
sensitive to angle of attack, with minimum 
damping occurring near zero degrees. 

For aeroelastic analysis of the wing, the 
mode shapes of the wing are required. The mode 
shapes were obtained from a finite element model 
which was based on measured wing stiffness and 
agreed well with wind-off node lines and 
frequencies measured during the tunnel test. 

The first eight mode shapes used in the 
aeroelastic analysis and their corresponding 
frequencies are shown in Figure 5. The first 
mode is first wing bending. The second mode Is 
primarily fore and aft motion with some second 
wing bending motion. The third mode is 
primarily second wing bending. The natural 
frequencies of the second and third modes are 
very close to one another. The fourth and fifth 
modes, first torsion and third bending 
respectively, also exhibit frequencies 
relatively close to one another. The same is 
also true of the sixth and seventh modes, which 
are second wing torsion and a coupling of second 
wing torsion and third wing bending. Further 
details of the wing and test program are given 
in ref. 11. 


XTRAN3S Program Description 

The modified unsteady transonic small 
perturbation (TSP) potential equation that is 
solved by XTRAN3S' is 

M 2 (<f> t + 2 »x>t = + F »x + G *y ] x 

v y x Y y ; y VY z ‘z (1) 

where the spatial coordinates x, y, and z are 
normalized by c r , the reference chord, and t 
is normalized by c r /V. Here the time scale 
factor k of ref. 7 is 1.0. The perturbation 
velocity potential $ is normalized by c r V. 
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The coefficients for the nonlinear terms in 


or 


2 U 


F = 

- \ (y+1)M 2 


G = 

\ (y-3)M 2 

(2) 

H = 

-(y-1)M 2 


F = 

- \ [3-(2-t)M 2 ]M 2 

(3) 

G = 

.Im 2 
2 " 


H = 

- M 2 



For F = G = H = 0 the linearized unsteady 
potential equation is obtained. The 
coefficients given by equation (2) are used for 
the calculations of this report for the ARW-2 
wing. 


The boundary conditions imposed on the 
outer edges of the computational region are 


♦ 

't’x + i’t 
♦z 


UzJ - Ux 


fy 

ftl 


= 0 upstream 

= 0 downstream 

= 0 above and below 

= 0 wing root 

= 0 far spanwise 

= 0 wake 


where [j indicates the jump in a quantity across 
the wake. 


The wing flow tangency condition is 


i = R" + 

r Z X 


z = 0" 


X L£ i. x i. 


TE 


where R = r/c r . 

The above conditions are those that are 
incorporated in the version of XTRAN3S code used 
for the computations reported in this paper 
(Version 1.5 with modified grid mapping). The 
outer boundary conditions recently have been 
improved by Whitlow 21 by implementing 
characteristic or " nonreflecting" boundary 
conditions. Calculations have indicated that 
these characteristic boundary conditions 
significantly reduce the reflection of 
disturbances from the computational boundaries. 

The XTRAN3S program has the capability for 
static and dynamic aeroelastic analyses of 
isolated wing configurations.' The structural 
equations of motion are based upon a linear 
analysis using natural vibration modes for both 
static and dynamic deformations. Numerical 
integration of the structural equations is 
performed for each time step along with the 
aerodynamic calculations giving the generalized 
coordinate of each natural mode at each time 
step. The natural vibration mode shapes and 
their associated generalized mass and 
generalized stiffness matrices are required as 
input data. 


Coordinate Transformation 


The finite difference grid contains 
60 x 20 x 40 points in the x, y, and z 
directions for a total of 48,000 points. In 
physical space the grid conforms to the wing 
planform and is mapped to a rectangular domain 
using the shearing transformation 

6 = ?(x,y), n = y, c = z, and t = t 


In computational space equation (1) becomes 

+ Fe xV + G5 y V + 2G Wn + G *n 2 


6 ( 4 ) 

* C7 Vs + *n } + VsVs + *n )] 

+ In Vs + V + %Vs + +n )] 

+ « t ^ 

The original version of XTRAN3S used 


S(x,y) = 


x - x LE (y) 


c(y) 


(5) 


to transform the streamwise coordinate. From 
equation (5), £ x = l/c(y), and Ey can be 
calculated analytically from equation (5). 

Using equation (5) results in a physical 
computational region whose streamwise extent is 
proportional to the local chord at each span 
station. For highly swept and tapered wings, 
the mesh is highly skewed in the far field of 
the physical domain. The resulting skewness of 
the grid led to numerical Instabilities that 
restricted the application of XTRAN3S to wings 
of low sweep and low taper. The grid has been 
revised such that in the physical domain, the 
computational region is a rectangular box thus 
alleviating the skewness of the grid. 

Versions of this type of grid arrangement 
have been given previously. 13 * 14 In ref. 13 a 
smooth stretching was used to map the regions 
from the wing to the upstream or downstream 
boundaries. Good results were obtained for the 
F-5 wing which is highly swept and highly 
tapered. In ref. 14, equation (5) is used on 
the wing and similar equations are used in the 
upstream and downstream regions with c(y) 
replaced by the local distance from the wing to 
the boundary of the region. This latter 
technique permits analytical evaluation of e x 
and 5y, but results in discontinuous values of 
these quantities at the leading and trailing 
edges. Even with these discontinuous values, 
results comparable to those of ref. 13 have been 
obtained for the F-5 wing. 14 

Herein, a finite difference grid is used 
that is similar to that of ref. 13. In the 
physical domain the grid is described 
analytically to give a smoothly-varying mesh 
spacing. The values of c x and ?y are then 
calculated numerically using second order finite 
difference formulae at each point in the 
computational domain. The c-distribution of 
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points (for all values of n) in the c-n domain 
is chosen to be the same as the x-dlstrlbution 
of points along the root chord in the physical 
domai n. 

Finite Difference Grid 


The x-y grid used for the computations 
presented in this report is shown in Figure 6 
along with the wing planform. The grid is 
defined in the following manner. First, 
equation (5) is used on the wing planform. The 
x-spacing of points along the chord essentially 
is that developed in references 12 and 22 , with 
39 points along the local chord; 38 points are 
equi spaced from x/c = 1/38 to 1.0, and an 
additional point is located at 
x/c = 1/(3 x 38). The grid extends from 20 wing 
root chords upstream of the wing root leading 
edge to 20 chords downstream of the wing root 
trailing edge. Eleven grid lines are used ahead 
of the wing and ten aft. The wing leading edge 
is centered between grid lines. The spacing 
ahead of the wing is given by 

x i u = X LE + A 1 ^ii + Vu ’ i u = 1 *” ,1 UP 

, X UP " X LE X 1 ‘ X LE 

3 = i (1 2 - 1 ) ' i 2 - 1 
1 UP ^ UP ' UP 

A 1 = X 1 " X LE ' A 3 

where 

xi = xj(y) = x-coordinate of first point 
on wing 


x LE = x LE (y) = x-coordinate of leading edge 


x MP = x-coordinate of upstream boundary 
(constant) 

and i u is the index of the points starting at 
the leading edge and running to the upstream 
boundary. Similarly downstream, 

x i . = X TE + Vd + B 3 1 d* 1 d = lf '** 1 DN 

a 

d - X DN " X TE X M ~ X TE 

3 1 (I 2 - 1) i 2 - 1 

1 DN' 1 DN ’DN 1 

B 1 = X W ' X TE ‘ B 3 

where 

x =x.,(y) = x-coordinate of first 
w " point downstream of wing 

x TC = x TC (y) = x-coordinate of trailing 
r£ IE edge 

x nM = x-coordinate of downstream boundary 
(constant) 


Fourteen rows of points are used along the 
span of the wing (including the points inboard 
of and on the plane of symmetry). The rows are 
distributed along the span using a cosine 
distribution 


- J + 1 


n j = n TIP 


c °s Cf ( w j . t -)]. J - 1. 


3 w 


where n-j-.p = semispan/c r and j w = number 
of grid rows on the wing (14 here). The first 
spanwise station is at 12 % span and the last one 
at 99.3% span. This distribution of grid points 
is used to emphasize the definition of lpads in 
the tip region which are critical for 
aeroelastic analyses. 


Outboard of the tip, the wing planform is 
extended smoothly to the far spanwise boundary 
at 1.5 semispan. The midchord line of the wing 
is extended to the far spanwise boundary as a 
parabola that has the slope of the wing midchord 
at the tip and is perpendicular to the outer 
grid boundary. The wing leading and trailing 
edges are extended using cubic equations that 
match leading and trailing edge slopes and 
intersect the outer grid boundary perpendicu- 
larly at one half the tip chord fore and aft of 
the midchord extension. The equations defining 
these grid points are 


x LE (n) = Xl E ti() + °i(n-n T ip) + D^n-n^) 
+ D 3 (n-n np ) 3 


where 


and 


dx 


°1 = 


LE 


I TIP 


An, 


1 2 C-ZDjAriQ + 3(X LE - X LE )] 


‘0 


OB 


‘TIP 


D 3 = r [D, Ao 0 - 2 (x, E - x. E )] 
Ang 1 U Lt 0B Lt TIP 




n 0B 


- n 


TIP 


n T jp = n-coordinate of wing-tip 


'OB 


= n-coordinate of far spanwise 
extent of grid 


x. _ = x-coordinate of wing-tip leading 

Lt TIP edge 


LE 


« x-coordinate of leading edge of 
OB outboard extension of wing 


Similarly for the trailing edge extension 

x te ( t i) = x TE TIp + E l^ n ' n TIP^ + E 2( n ' n TIP^ 
+ E3(n-n TIp ) 


and id is the index of the points starting 
aft of the trailing edge and running to the 
downstream boundary. 
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where 


E, 

= ^TE 


1 

dn 

TIP 

c 

1 

t-2E«Aru + 3(Xjr * Xyr )3 
1 v fc 0B TIP 

E 2 

i 

> 

O 

c 

= 1 

[£i Ann • 2 ( x tc ~ X TF ^ 

i v it QB it TIp 

E 3 



and 

x TF = x-coordinate of wing trailing 
b TIP edge at tip 

x T p = x-coordinate of trailing edge of 
^OB outboard extension of wing 

The outboard extension of the midchord line 
intercepts the right hand boundary at 


= 


"OB 


TIP 


, dx M 

l dn 


TIp) An 0 


where 



= (Dj + E x )/2 


and offsetting by ±c T j p /2, 


XLE 0B = XM 08 ' CTIp/2 
XTE 0B = \)B + ° TIp/2 


The same streamwise distribution of points that 
is used on the wing is used between these 
outboard extension of the leading and trailing 
edge grid lines. 

The n-distributlon of points outboard of 
the wing-tip is given by a cubic equation. The 
first row of points outboard of the wing and the 
last row on the wing are located symmetrically 
about the tip. Also applying a zero second 
derivative condition at the tip gives 

n jo = C lJo + C 3 j o + n TIP ’ J o = 1 * j 0B 

r _ j 0B (n 01 ~ n TIP^ ' (n 0B ~ n TIP^ 

- ’ 2 

dg B (j 0B " 1) 

c _ ^ n 0B " n TIp) ' jpB^Ol ' ^TIP^ 
jQB^OB * ^ 


= number of spanwise points outboard of 
UB wing-tip 

j = index of points outboard of wing-tip 

n m = n-location of first point outboard of 
u wing-tip 

In the z-direction, the grid extends 25 
root chords above and below the wing. Twenty 
planes of points are used above and 20 planes 
are used symmetrically below the wing. The 
distribution of points is that developed and 
applied in references 12 and 22. 

The resulting grid Is smooth and embeds the 
planform smoothly. This grid appears to be 
reasonable, but since to date only limited 
variations have been Investigated for three 
dimensional configurations, the grid cannot be 
considered an optimal choice. The grid is . 
reasonably fine on the wing but is relatively 
coarse off the wing. 


R e suits and Discussion 


Calculations have been made with XTRAN3S 
for the DAST ARW-2 wing for a range of Mach 
numbers from 0.60 to 0.90, for zero wing root 
angle of attack, and for the aileron 
undeflected. Static aeroelastic calculations 
for the flexilbe wing have been made for dynamic 
pressures up to 1000 psf. Static pressure 
calculations were made for the rigid wing for a 
limited range of Mach numbers. Limited dynamic 
aeroelastic calculations were made to explore 
the wing instability near M = 0.90. 

The measured airfoil ordinates have been 
used to determine the wing surface slopes 
required for input to XTRAN3S. The measured 
ordinates were fitted with a parametric spline 
with smoothing , zi and the slopes were calculated 
from the spline fit. The mode shapes from the 
finite element model were Interpolated to the 
XTRAN3S computational grid locations with a 
surface spline. 24 

Steady Pressure Distributions .- For the rigid 
wing, calculated static pressures for several 
Mach numbers are shown in Figure 7. The upper 
and lower surface pressures are shown for zero 
angle of attack. The upper surface pressures 
are supercritical at M = 0.75 in the outboard 
region of the wing, and at M = 0.80 a strong 
shock wave is calculated which is very 
pronounced In the outboard region of the wing. 
Attempts to obtain the static rigid solution at 
M = 0.81 were unsuccessful. 

The variation of calculated static 
pressures with dynamic pressure at M = 0.8 for 
the flexible wing is shown In Figure 8. The 
dynamic pressure is varied from zero (rigid 
wing) to 900 psf. The addition of wing 
flexibility alleviates the strong shock wave on 
the upper surface in the outboard region of the 
wing. This is a result of the wing-tip washout 
which occurs with increasing dynamic pressure 
wherein the wing is twisted nosedown and 
deflected upwards. 
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Calculated static pressures for the 
flexible wing for several Mach numbers and a 
dynamic pressure of 300 psf are shown in Figure 
9. With the addition of the static aeroelastic 
deformation, calculations with XTRAN3S were made 
up to M = 0.90. The upper surface pressures are 
supercritical at M = 0.80 with the shock being 
the strongest in the center span region. With 
Increasing Mach number the upper surface shock 
becomes stronger and moves to the trailing 
edge. In addition, a shock begins to form on 
the lower surface in the wing-tip region at 
M = 0.85 and becomes stronger and moves Inboard 
at M = 0.90. The strong shocks located aft on 
the wing for M = 0.85 and 0.90 are not 
considered realistic and these cases are 
considerably beyond the design Mach number of 
0.80. The pressure distributions at these 
conditions are similar to those seen In 
two-dimensional potential flow calculations when 
the assumption of isentropic flow is 

violated. 5 

Comparisons of calculated static pressures 
for the flexible wing with wind tunnel data are 
shown in Figure 10 for M = 0.70, 0.80, and 
0.85. Figure 10(a) shows upper and lower 
surface static pressure distributions for 
M = 0.70, a dynamic pressure of 200 psf and zero 
angle of attack. The agreement of the 
calculated pressures with experiment is only 
fair. The lower surface pressure in the Inboard 
region is underpredicted along the chord and 
there is a consistent underprediction in the 
leading edge region along the span. The 
agreement between theory and experiment is 
better in the outer span region. Although not 
Shown here, similar agreement between measured 
and calculated pressures was observed for 
q = 100 and 300 psf at M = 0.70 and for q = 100, 
200, and 300 psf at M = 0.60. The flow is 
subcritical at these conditions. 

Similar agreement is not seen for cases 
with supercritical flow. Figures 10b, 10c, and 
lOd show comparisons of measured and calculated 
pressures at M = 0.8 for q = 100, 200 and 300 
psf, respectively. No shocks are seen in the 
measured data for these cases. In contrast, the 
calculations show moderate to strong shocks 
except for the outboard station at q = 200 
(Fig. 10c) and 300 psf (Fig. lOd). Also, the 
calculations overpredict the upper surface 
pressures and underpredict the lower surface 
pressures on the forward portion of the wing. 
Increasing dynamic pressure increases wing-tip 
washout and the agreement with experiment 
improves in the wing-tip region. For q = 200 
psf, (Figure 10(c)), the shock calculated at the 
outboard station is in good agreement with 
experiment. As the dynamic pressure Increases 
to 300 psf, (Figure 10(d)), the shock vanishes 
outboard. Finally, Figure 10(e) shows the 
static pressure distributions at M = 0.85 for 
q = 200 psf. Similar to M = 0.8, the upper 
surface shock strength is overpredicted over the 
entire wing and the lower surface pressure is 
underpredicted inboard and along the wing 
leading edge. 

The calculations using XTRAN3S have shown 
reasonable trends but only fair agreement with 
the measured data for cases of subcritical 
flow. Larger discrepancies are indicated for 


cases with strong calculated shock waves. The 
calculations do not include either the effects 
of the fuselage or a viscous boundary layer, 
each of which may have a large effect in 
transonic flow. Further effort is needed to 
improve the agreement with the experimental data 
possibly by including the fuselage and boundary 
layer, or by improving the analysis in other 
respects such as the number of grid points and 
so forth. 

Static Aeroelastic Effects. - The static 
aeroelastic deformation at the wing-tip 
is shown in Figure 11. The experimental and 
calculated wing-tip vertical displacement, zt, 
and rotation, e^, are shown for varying 
dynamic pressure at M = 0.8 and zero wing root 
angle of attack. The calculations were made for 
increments of 100 psf in q. Roth the 
computational and experimental wing-tip 
deflections, zt, are nonlinear with dynamic 
pressure which would not be expected from 
linearized aerodynamic theories. The wing-tip 
displacement is approximately 40% overpredicted 
at q = 300 psf, which is expected from the 
overpredicted upper surface pressures shown in 
Figure 10(d). The calculated tip displacement 
reaches a maximum of approximately 10 inches at 
q = 600 psf. The calculated wing tip twist 
angle, 9^, is in good agreement with the 
experimental data although this agreement may be 
fortuitous. One test point taken in air is 
included in Figure 11 and agrees with the 
measurements in Freon within the range of 
experimental scatter. 

A comparison of calculated wing-tip 
deflections at two Mach numbers versus dynamic 
pressure is shown in Figure 12. The predicted 
deflection data at M = 0.70 are very similar to 
that for M = 0.80. At the lower Mach number the 
wing-tip deflections are also nonlinear with 
dynamic pressure even though the flow is 
subcritical for M = 0.70. Although the tip 
deflections are slightly less at M = 0.70, they 
are of the same order of magnitude as at 
M = 0.80. 

The variation of calculated wing-tip 
deflections with Mach number at q = 300 psf is 
shown in Figure 13. The calculations were made 
for increments Of 0.05 in Mach number up 
M = 0.80 and for increments of 0.01 between 
M = 0.80 and M = 0.90. The wing-tip vertical 
displacement and rotation vary slightly with 
Mach number between M = 0.60 to 0.80. However, 
above M = 0.80 the wing-tip nose-down twist 
increases rapidly with increasing Mach number 
and the wing displacement rises sharply to a 
maximum. These effects are caused by the 
development of strong shocks moving aft along 
the wing as seen for M = 0.85 and 0.90 in Figure 
9. The sharply increasing nose down twist of 
the wing-tip with increasing Mach number above 
M = 0.80 is qualitatively similar to that 
observed on the model during the test. 

Similarly, the wing-tip displacement shown in 
Figure 13 is similar to that observed in the 
test even though the detailed calculated 
pressures do not compare well with experiment. 
The rapid Increase in wing-tip displacement and 
rotation beyond M = 0.80 is expected for Mach 
numbers above the design cruise condition. 
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The variation of predicted wing lift and 
pitching moment with Mach number for the 
flexible (q = 300) and rigid wing (q = 0) is 
shown in Figure 14. The results for both the 
rigid and flexible wing show a nonlinear 
increase in the magnitude of the lift and the 
moment coefficients with Mach number. As 
previously noted, attempts to calculate the 
rigid wing above M = 0.80 were unsuccessful. 
Above M = 0.80 the forces on the flexible wing 
rapidly increase until M « 0.88 where the loads 
reach a maximum and begin to decrease. 

The variation of predicted wing lift and 
pitching moment with dynamic pressure for 
M = 0.70 and M = 0.80 is shown in Figure 15. As 
in the comparison of wing-tip deflections at 
these Mach numbers, the data for M = 0.70 is 
very similar to that for M = 0.80. For both 
Mach numbers the effect of wing flexibility is 
quite significant. 

The trends shown here for the effects of 
flexibility are large, significant, and are of 
considerable interest. It nust be kept in mind 
that these results will be affected by the lack 
of agreement of the calculated and measured 
static pressures. Further effort is needed to 
calculate these efforts accurately with codes 
such as XTRAN3S on thick supercritical wings. 

Aeroelastic Stability Calculations .- Dynamic 
aeroelastic calculations were made with XTRAN3S 
in an attempt to predict the instability 
boundary encountered during the wind tunnel test 
and shown in Figure 4. The aeroelastic 
transient calculations were made for q = 300 psf 
at Mach numbers up to 0.90. No instability was 
encountered. It should be noted that, as 
previously shown, the agreement of calculated 
pressure distributions with experiment degrade 
rapidly for Mach numbers greater than 0.70. If 
the angle of attack were adjusted to match 
pressures at the higher Mach numbers, it is 
possible that the instability boundary might be 
calculated. 


Concluding Remarks 

An overview of the transonic small 
perturbation program XTRAN3S and a 
recently-developed grid coordinate 
transformation has been given. Results 
calculated using XTRAN3S have been compared with 
data measured for the DAST ARW-2 wing in the 
Langley Transonic Dynamics Tunnel. Several 
static aerodynamic cases were calculated for a 
rigid wing and with flexibility included. The 
addition of wing flexibility was shown to 
strongly influence the predicted static pressure 
for this wing. Predicted static pressures for 
the flexible wing were compared to wind tunnel 
data for Mach numbers from 0.70 to 0.85 for 
dynamic pressures from 100 to 300 psf and for 
zero wing root angle of attack. The calculated 
pressures were in fair agreement with 
experimental data for subcritical flow cases. 

For cases where the flow was supercritical the 
program overpredicted transonic effects such as 
the shock wave strength. The pressure along the 
leading edge was poorly predicted. The best 
agreement of pressures occurred near the 
wing-tip. 


Trends of predicted wing-tip deflection and 

wing forces for a flexible wing were studied. 

At M = 0.80 the predicted wing-tip displacement 
was approximately 40% high, but the tip rotation 
was close to the experimental value. The 
program predicted trends for wing-tip deflection 
which were qualitatively similar to those 
observed during the wind tunnel test, but were 
quantitatively overpredicted. The wing-tip 
deflections and forces were shown to be 
nonlinear with dynamic pressure and Mach 
number. The predicted displacements and forces 
rapidly increased in magnitude beyond the design 
Mach number of 0.80 and peaked at approximately 
M = 0.88. 

The XTRAN3S program was also used to 
investigate the dynamic stability of the wing. 

An attempt was made to calculate the instability 
observed during the tunnel test at M « 0.90. 

This attempt was unsuccessful and further effort 
is needed to improve the accuracy of the static 
flow field predictions which may in turn improve 
the calculation of this instability. 
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Fig. 1 Wing mounted In TDT test section. 



Mach number 


Fig. 3 Wind tunnel test conditions. 




.5 .6 .7 .8 .9 i.o 


Mach number 


Fig. 2 Wing planform and instrumentation 
locations <1n inches). 


Fig. 4 Measured wing instability boundary. 
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b) (continued) M = 0.80, q = 100 psf. 


Fig. 10 Comparison of measured and calculated steady pressure distributions 
for the flexible wing, a = 0°. 
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c) (continued) M = 0.80, q = 200 psf. 
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e) (concluded) M = 0.85, q = 200 psf. 
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Fig. 11 Comparison of measured and calculated 
wing-tip displacement and twist versus 
dynamic pressure, M = 0.8, a = 0°. 


Fig. 12 Calculated wing-tip displacement and 
twist versus dynamic pressure, a = 0 
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Fig. 14 Calculated wing lift and pitching 

moment coefficients versus Mach number, 

a = 0°. 


Fig. 15 Calculated wing lift and pitching 
moment coefficients versus dynamic 
pressure, a = 0°. 
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